Fano blockade by a Bose-Einstein condensate in an optical lattice 
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We study the transport of atoms across a localized Bose-Einstein condensate in a one-dimensional 
optical lattice. For atoms scattering off the condensate we predict total reflection as well as full 
transmission for certain parameter values on the basis of an exactly solvable model. The findings 
of analytical and numerical calculations are interpreted by a tunable Fano-like resonance and may 
lead to interesting applications for blocking and filtering atom beams. 
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An understanding of the transport properties of ultra- 
cold atoms is vital for the development of technologi- 
cal applications in the fields of matter-wave interferom- 
etry [l| or quantum information processing with neutral 
atoms 0, H 0, H| • Over the last couple of years it has 
been shown that optical lattices, generated by counter- 
propagating laser beams and providing a periodic poten- 
tial modulation for the atoms, introduce many interesting 
and potentially useful effects by modifying single atom 
properties and enhancing correlations between atoms . 
Here we discuss the scattering of atoms across a localized 
Bose-Einstein condensate (BEC) in an optical lattice. We 
find that dramatic effects of scattering resonances with 
either full transparency or total reflection can occur. 

Previously, transparency effects have been conjectured 
for the scattering of He atoms on a film of superfluid 
helium-4 and similar effects have been predicted for 
the scattering of atoms on a BEC in a trap of finite depth 
These effects were first attributed to the coherent 
interactions within the target and with the scattering 
atoms but a full understanding of the numerical results 
was not achieved. More numerical results were produced 
later and a Levinson theorem was proved on general 
grounds without revealing the mechanism for trans- 
parency effects. 

In this Letter we present and analyse a very simple 
and analytically solvable one-dimensional model of atom 
scattering by a BEC. The model shows transparency as 
well as blockade of atoms by total reflection, which is 
interpreted as a Fano resonance. In the problem stud- 
ied by Fano, the resonance may both enhance and sup- 
press scattering due to interference Here the atom- 
atom interaction leads to an effective nonlinearity. It was 
shown recently |l2j. that nonlinearity generates several 
scattering channels which can lead to destructive inter- 
ference and Fano resonances similar to the original Fano 
problem. Proposed applications in nonlinear optics and 
Josephson junction networks encounter various difficul- 
ties due to inhomogeneities and dissipation 0|. They 
are absent in the present study, thus making the reso- 
nant atom-BEC scattering ideal for harvesting on Fano 
resonances. 



FIG. 1: (Color online). Scattering scheme in an optical 
lattice. The incoming, reflected, and transmitted beams of 
atoms are represented as plane waves. The BEC is centered 
in n — n c , the only nonlinear lattice site. 



We consider a BEC on a lattice, where interactions 
between atoms are present on one lattice site only (see 
Fig. [Q. Such a situation could be realized experimen- 
tally by combining optical lattices with atom-chip tech- 
nology 0, 0| or in optical micro- lense arrays |l6j where 
the s-wave scattering length of atoms can be tuned 
by an inhomogeneous magnetic 0, or laser field 
|20| . Specifically we consider the discrete nonlinear 
Schrodinger (DNLS) equation, a classical variant of the 
Bose-Hubbard model appropriate for a BEC in a periodic 
potential in the tight binding limit 6] . With interactions 
being present only on site number n c , we write in dimen- 
sionless form 

= - + - Tl^rJ^nA.^, (i) 

at 

where $„(t) is a complex amplitude of the BEC field 
at site n and —7 = U/J is the interaction strength on 
site n c . This simple model reflects generic features of 
BECs in a one-dimensional optical lattice with inhomo- 
geneous scattering length. Furthermore, this model could 
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be realized quantitatively in a deep optical lattice with 
tight transverse confinement pi) . For atoms with mass 
M in a lattice with spacing d in the tight binding limit, 
J sa 4s 3 / 4 e~ 2v/ ^ / y/nE r is the energy scale for tunneling 
between the lattice sites, where s = Vo/E r is the depth 
of the optical lattice Vb measured in units of the recoil 
energy E r — 2fm 2 j (d 2 M) . The on-site interaction en- 
ergy per atom is U = 4Tra s h 2 J d 3 x\ip(x)\ 4: /M , where a s 
is the tunable s-wave scattering length at the nonlinear 
site n c and ^(x) is the localized Wannier state associated 
with the lowest Bloch band of the lattice. The number 
of atoms in the lattice is given by N = l^n| 2 - Small- 
amplitude plane-wave solutions of Eq. Q take the form 
= exp(ikn) exp(-iEkt) and satisfy the relation 

£ fc = -2cosfc, (2) 

which defines the band of single-particle energies Ek S 
[—2,2] [see FigJ^a)]. The unit of dimensionless energy 
is J and the quasi-momentum k is measured in units of 

d-\ 

First, we look for localized and stationary solutions of 
Eq.JTJ, corresponding to the BEC centered in n = n c . 
For simplicity we assume that interactions are attractive 
and thus 7 > 0. As an ansatz, we take an exponentially 
localized profile: = b n (t) = b nc x^ n ~ n ^ exp(—iE b t), 

where b nr is the condensate amplitude, |x| < 1, and E b is 
the respective energy. Inserting this expression into ijTJl. 
we obtain that 

E b = + and x = -(E b + g)/2, (3) 

where g = -fb 2 lc (g > 0). Eqs. © correspond to solu- 
tions for localized BECs with E b being outside of the 
band Ek [Eb < —2, see Fig^a)]. They are similar 
to bright lattice solitons pinned to the nonlinear lattice 
site. These localized modes exist only above a threshold 
N b = -E b /-y > Nl h H3, given by N* h = 2/7. We as- 
sume that Nb is significantly larger than the threshold, 
which should be easily achieved in a possible experiment. 

By controlling the number of atoms in the BEC or by 
tuning the nonlinear coefficient, we can easily modify the 
BEC energy E b = —jN b (this is equivalent to modifying 
the parameter g). This is one of the keys for a tunable 
Fano-blockade scheme. As we will show later, this energy 
is directly related to the energy where zero transmission 
of the atom beam through the BEC is observed. 

We consider three different values for g as examples: 
gi = 0.36, <72 = 0.6, and 53 = 0.9. The corresponding 
energies [E b \ = -2.03 (box), E b 2 = —2.09 (diamond), 
and Ebs = —2.19 (triangle)] and the profiles of the BEC 
states are shown in the left part of Fig. |5fb). By do- 
ing a linear stability analysis, we find that all these so- 
lutions centered in n — n c are stable. This is impor- 
tant for two reasons. First, experimental generation of 
these states should be possible. Second, the scattering of 
small-amplitude beams by these modes can be viewed as 
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FIG. 2: (Color online), (a) Energy diagram for localized and 
extended solutions. In the linear band, Ek is plotted, (b) 
Zoom of the region E ~ —2. Boxes, diamonds, and triangles 
correspond to the BEC (Eb < —2) and to the local mode 
(El > —2) solutions for gi, Q2, and gz, respectively. The 
corresponding BEC profiles are shown. The filled circles cor- 
respond to the resonance condition 111 2t . 

a small perturbation of the latter, which will stay small if 
the BEC is stable, but would grow (and the BEC would 
decay) otherwise. 

In order to study the scattering of a propagating atom 
beam by the localized BEC, we consider: 

^n(t) = K{t)+bn{t), (4) 

where <fi n -C b n . We linearize Eq.lJIJ with respect to 
4> n {t), obtaining the following equation for the atom 
beam: 

(5) 

Far away from n — n c , the solution of this equation 
corresponds to a propagating plane wave which satisfies 
(J2J. The localized BEC generates a scattering potential 
for the atom beam which has a constant and a time- 
dependent part. Equation (Jjjjl is similar to the one found 
111 Ref.0 for the scattering of plane waves against a dis- 
crete breather. In that work the appearance of Fano res- 
onances was studied in a full nonlinear lattice in the 
approximation of a very localized nonlinear mode. The 
scattering potential is then reduced to one lattice site and 
implies j> 1, which makes the discrete breather a very 
strong scatterer for almost all incoming plane waves. Our 
setup allows us to tune g to smaller values, and thus ad- 
mits Fano resonances on a background of almost perfect 
transmission. 

We solve the equation (JSJ by using a Bogoliubov trans- 
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formation Q given by 



-iEt , * -i(2E b -E)t 



(6) 



and by inserting it in JSJ) we obtain the discrete Bogoli- 
ubov (DB) equations: 

Eu n = -(u n+ i + u n _i) - g(2u na + v„ c )S n<nc , (7) 
(2E b - E)v n = -(v n+ i + v n -i) - g(2v na + u„J^,„ c .(8) 

Here u n corresponds to an open channel which, far away 
from n C) represents a propagating atom beam for which 
the energy is in the band E = E^ £ [—2,2]. Contrary, 
v n represents a closed channel whose extended states far 
away from the scattering center have 2Eb — E ^ [—2, 2]. 
They are thus located outside the open channel contin- 
uum and cannot be excited in the same energy range. 
However, even in the absence of any coupling between 
both channels, the scattering center provides also a lo- 
calized state in the spectrum of the closed channel with 
energy El- The localized state, by definition, is located 
outside the band of extended states in the v n channel, but 
may be located inside the u„-channel band E £ [—2,2]. 
In such a case, taking the coupling between channels into 
account, we encounter a Fano resonance for El = E^. 

Let us consider the situation when both channels are 
decoupled. For this particular situation, the closed chan- 
nel equation is given by 



■(v n +i + v n -i) - 2gv n J n 



and admits a localized solution: v„ = v„ w' n ™ c ' (\w\ < 



(2E b - E)v n 

idmits a 
1), for which 



w = — g 



and E = El = 2(E b + ^1+ g 2 ). 



We call this solution the local mode (LM). El corre- 
sponds to the LM energy, which is always inside the con- 
tinuum of the open channel: if g — > El — > —2 and, 
if g — > oo El — > 0. Therefore, due to the time depen- 
dence of the original scattering potential, the closed chan- 
nel is able to resonate with the open one at a frequency 
that depends on externally controllable parameters. 

Keeping in mind the localized nature of the LM and 
the propagating one of the open channel, we make the 
following ansatz: 



oi e lk(n - n "'> + bi e - 2fe ("-™c) ; n < n c 

lk(n-n c ) . n> n ) I 9 ) 



ci e 



\n—n c \ 



(10) 



Here, a\, bi and ci, represent the incoming, reflected 
and transmitted beam amplitudes, respectively. v nc cor- 
respond to the closed channel amplitude and \w\ < 1. 
The beam quasimomcntum k can be generated in the ex- 
periment by using a phase imprinting method [23 | , Bragg 
scattering, or simply by acceleration of the matter- wave 
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FIG. 3: T versus k and g. The dotted white line correspond 
to T(fci) = 0, the resonance position. 



probe in an external potential. We solve analytically the 
scattering problem by inserting © and (|10|l in (JJJ) and 
(JHJ for n — n Cl n c ± 1. By doing so, we obtain that the 
open channel satisfies P|l. en + 61 = ci . w = w. and that 
the transmission is given by: 



T{k) 



4 sin 2 k 



4 sin 2 k + ( 2g 



^/(E k -2E b ) 2 -A-2g 



(11) 



(T = |ci/ai| 2 ). Resonances occur when the denomina- 
tor diverges or when yJ(E k - 2E b ) 2 - 4 - 2g = 0. The 



condition for the resonance is 



Ek = E L => k L = arccos(— El/2). 



(12) 



Ea. (|12|) implies that the transmission for atom beams 
through the BEC is reduced to zero when a LM is gen- 
erated in the process, i.e. when the closed channel res- 
onates with the open one. 

In Fig[5Jb), we show the corresponding energies for 
the LM's, for the three values of g, Eli — —1-94 (box), 
El2 = —1-84 (diamond), and El3 — —1.69 (triangle). 
The curve in that figure corresponds to Eq. (0), and the 
filled circles correspond to the value of k for which the 
resonance is expected ()12|l . In principle, any energy in 
the interval {—2, 0} can be a good candidate for the ob- 
servation of a Fano resonance in this setup. But, as we 
will show later, the response of the system is not always 
the same and it essentially depends on the BEC profile. 

We compute the transmission T = T(k, g) and we vary 
the beam velocity k and the nonlinear parameter g. FigEl 
shows this sweep of parameters, where dark and bright re- 
gions represent a low (T — ► 0) and a high (T — > 1) trans- 
mission response, respectively. In Fig^a) we show some 
T(k) for three values of g. As g increases the width and 
the position (dashed white line in Fig.^J of the resonance 
increase. Thus, the more localized the BEC becomes, the 
stronger it reflects the atom beam off resonance. As ex- 
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three sites of the array, similar transmission curves are 
observed. This shows the robustness of our theoretical 
results in a more realistic experimental setup. 

In conclusion, we have investigated Fano resonances 
in the context of Bose-Einstein condensates in an optical 
lattice. The implementation of this idea can be viewed as 
a powerful tool for controlling the transmission of matter 
waves in interferometry and quantum information pro- 
cesses. Fano resonances rely on destructive interference 
and are thus inherent to wave dynamics. An observa- 
tion of these resonances in atom-BEC scattering would 
provide, in addition to tunable filters, a new demonstra- 
tion of the quantum matter wave character of ultracold 
atoms. 



60 120 150 300 60120 
t t t 



FIG. 4: (Color online), (a) T versus k for the three values of 
g. Lines: eq. Hill , points: real time numerical simulations of 
eq. Jy for gi (boxes), gi (diamonds), and (73 (triangles), (b)- 
(d) Evolution of |^ n (t)| 2 in space and time: (b) gi,k = 1.37, 
(c) 32, A; = 0.39, (d) g 3y k = 1.37. 

pected from (|12fl , increasing g (which corresponds to de- 
creasing Eb or increasing Nt) leads to an increase of the 
resonance energy. By tuning the nonlinear parameter g, 
we can thus choose the amount of the beam which passes 
through the BEC. Off resonance (for larger values of k), 
we can select the percentage of the incoming beam that 
is transmitted for a defined quasimomcntum. Therefore, 
the actual setup can be used as a 100% blockade or as a 
selective filter. 

Now, we look for a numerical confirmation of our the- 
oretical description for the scattering of an atom beam 
against a localized BEC. In the simulations, we ini- 
tialize the atom beam with a Gaussian profile: <j) = 
4>q exp[— a(n — n ) 2 } exp[i/c(n — no)] with (j) Q — 0.01 and 
a — 0.001. n is the initial location of the center of the 
distribution and well separated from the BEC to avoid 
an initial interaction, k is the initial quasimomcntum. 
The amplitude 4>q was chosen to be very small compared 
to the BEC amplitude {4>o/b nc ~ 1%) in order to justify 
©. The value of a implies a spatial width of approxi- 
mately 60 sites and a reciprocal width in fc-space of 0.12. 
With this choice we can clearly observe the resonant re- 
sponse of the system. In Fig^Ja) the symbols denote 
our numerical results for three different values of g. The 
agreement between theory and simulations is almost per- 
fect. We have some disagreement for small values of k 
where the group velocity is too small and the computa- 
tion of T is more complicated. FigsQ](b), (c), and (d) 
show some numerical examples of the scattering process 
with a transmission of 80%, 0% and 40%, respectively. It 
is worth mentioning that, if we extend the nonlinearity to 



[1] T. Schumm et ai, Nat. Phys. 1, 57 (2005). 

[2] D. Jaksch et al, Phys. Rev. Lett. 82, 1975 (1999). 

[3] G. K. Brennen, C. M. Caves, P. S. Jessen, and I. H. 

Deutsch, Phys. Rev. Lett. 82, 1060 (1999). 
[4] J. K. Pachos and P. L. Knight, Phys. Rev. Lett. 91, 

107902 (2003). 

[5] A. Kay, J. K. Pachos, and C. S. Adams, Phys. Rev. A 

73, 022310 (2006). 
[6] O. Morsch and M. Oberthaler, Rev. Mod. Phys. 78, 179 

(2006). 

[7] J. W. Halley, C. E. Campbell, C. F. Giese, and K. Goetz, 

Phys. Rev. Lett. 71, 2429 (1993). 
[8] A. Wynveen et ai, Phys. Rev. A 62, 023602 (2000). 
[9] U. V. Poulsen and K. M0lmer, Phys. Rev. A 67, 013610 

(2003). 

[10] J. Brand, I. Haring, and J.-M. Rost, Phys. Rev. Lett. 
91, 070403 (2003); T. Cretegny, S. Aubry, and S. Flach, 
Physica D 119, 73 (1998). 

[11] U. Fano, Phys. Rev. 124, 1866 (1961). 

[12] S. Flach, A. E. Miroshnichenko, V. Fleurov, and M. V. 
Fistul, Phys. Rev. Lett. 90, 084101 (2003); A. E. Mirosh- 
nichenko, S. Flach, and B. Malomed, CHAOS 13, 874 
(2003). 

[13] S. Flach, V. Fleurov, A. Gorbach, and A. E. Mirosh- 
nichenko, Phys. Rev. Lett. 95, 023901 (2005); A. E. 
Miroshnichenko et ai, Phys. Rev. B 71, 174306 (2005). 

[14] W. Hansel, P. Hommelhoff, T. W. Hansch, and J. Reichel, 
Nature 413, 498 (2001). 

[15] H. Ott et al, Phys. Rev. Lett. 87, 230401 (2001). 

[16] R. Dumke et al, Phys. Rev. Lett. 89, 220402 (2002). 

[17] E. Tiesinga, B. J. Verhaar, and H. T. C. Stoof, Phys. 
Rev. A 47, 4114 (1993). 

[18] S. Inouye et ai, Nature 392, 151 (1998). 

[19] P. O. Fedichev, Y. Kagan, G. V. Shlyapnikov, and 
J. T. M. Walraven, Phys. Rev. Lett. 77, 2913 (1996). 

[20] M. Theis et ai, Phys. Rev. Lett. 93, 123001 (2004). 

[21] A. Smerzi and A. Trombettoni, Phys. Rev. A 68, 023613 
(2003). 

[22] G. P. Tsironis, M. I. Molina, and D. Hennig, Phys. Rev. 

E 50, 2365 (1994). 
[23] J. Denschlag et al, Science 287, 97 (2000). 



